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Abstract The centrality dependence of transverse momentum spectra for identified hadrons at midrapidity 
in Au+Au collisions at ^/sjvjv = 200 GeV is systematically studied in a quark combination model. The pt 
spectra of 7T , K , p(p) and A(A) in different centrality bins and the nuclear modification factors (Rcp) for 
these hadrons are calculated. The centrality dependence of the average collective transverse velocity (P(r)) for 
the hot and dense quark matter is obtained in Au+Au collisions, and it is applied to a relative smaller Cu+Cu 
collision system. The centrality dependence of pt spectra and the Rcp for 7r°, and A in Cu+Cu collisions 
at tJsnn = 200 GeV are well described. The results show that {/3 (r)) is only a function of the number of 
participants N vart and it is independent of the collision system. 
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1 Introduction 

The hadronic collectivity is one of the important 
properties in ultra-rclativistic heavy ion collisions. It 
provides a lot of information on the initial spatial 
anisotropy of the reaction zone^'— '— '^ , the degree of 
thermalization ^'^ and the hadronization mechanism 
of the hot and dense medium produced in nucleus- 
nucleus collisions' 2 ^. It can also help to understand 
the broadening of jetlike particle correlations^ and 
high p T jet-quenching^'^. Strong partonic mul- 
tiple scatterings in nucleus-nucleus collisions would 
generate the collectivity of quarks, which then re- 
sult in the observed collectivity of final hadrons. The 
quark number scaling of the hadronic elliptic flow u 2 
is a piece of evidence for this original quark collec- 
tivity li£'2&iiJ. Being a key hadronization mech- 
anism, quark combination picture has successfully 
described many features of multi-particle produc- 
tion in high energy heavy ion collisions, e.g. the 
high p/ir ratio in intermediate transverse momentum 
region^'— '^, the quark number scaling behavior of 
hadron elliptic flow^— '^J and its fine structure at 



small Pt^j hadron longitudinal and transverse mo- 
mentum distributions' 2 '—'—'^, the yields and mul- 
tiplicity ratios^. Therefore, one can extract the col- 
lectivity of the hot and dense quark matter from the 
data of hadrons through quark combination mecha- 
nism. 

In general, the nucleus-nucleus collisions in dif- 
ferent centralities will produce different sizes of the 
hot and dense quark matter. This would cause the 
transverse collective flow for the hot and dense quark 
matter varying with collision centralities. This vari- 
ance of collective flow in quark level would be em- 
bodied in the transverse momentum spectra of fi- 
nal hadrons'^' ^ and particularly in their nuclear 
modification factors Rcp- The thermal and hydro- 
dynamic models have described the centrality de- 
pendence of the transverse momentum distributions 
for final hadrons in low p T region with a statisti- 
cal hadronization method' 22 ' — ' — ' — ' — ' ^ . Using the 
quark recombination at intermediate p T and parton 
fragmentation at high p T , the Duke group in Ref. 
™ has explained the strong suppression of hadron 
Rcp at high transverse momenta, and the baryon- 
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meson difference of hadron R C p in intermediate p T 
region, fn the present paper, with quark combina- 
tion at all p T , we use our quark combination model 
to study the hadron p T spectra from central to pe- 
ripheral collisions. We investigate the p T spectra of 
identified hadrons at midrapidity in different central- 
izes in Au+Au collisions at ^Js NN = 200 GeV to ob- 
tain the centrality dependence of transverse collective 
flow for the hot and dense quark matter. Further- 
more, we apply it to relative smaller Cu+Cu collision 
systems, calculate the transverse momentum spectra 
of final hadrons and nuclear modification factors Rcp 
in Cu+Cu collisions at ^Js NN = 200 GeV, and com- 
pare them with the experimental data from STAR 
and PHEN1X Collaborations. 

In the next two sections, we briefly introduce the 
quark combination model and the transverse momen- 
tum spectra of quarks just before hadronization. The 
results and discussions are in Sect. 4. Summary is 
given in Sect. 5. 



2 The quark combination model 



Within the same quark combination mechanism, 
all kinds of combination-like models, such as recom- 
bination model^'— '^J, and coalescence model li&iSJj 
have their own features. Our quark combination 
model was first proposed for high energy e + e~ and 
pp collisions'^'—'—'—'—'^'. Recently we have ex- 
tended the model to ultra-relativistic heavy ion 
collisions^'— '—'^J. The model describes the pro- 
duction of initially produced ground state mesons 
(36— plets) and baryons (56 — plets). In principle, it 
can also be applied to the production of excited states 
^ and exotic states^ . These hadrons through com- 
bination of constituent quarks are then allowed to 
decay into the final state hadrons. We take into ac- 
count the decay contributions of all resonances of 
56 — plet baryons and 36 — plet mesons, and cover 
all available decay channels by using the decay pro- 
gram of PYTHIA 6.1 The main idea is to line up 
quarks and anti-quarks in a one-dimensional order in 
phase space, e.g. in rapidity, and let them combine 
into initial hadrons one by one following a combina- 
tion rule. See the second section of Ref. ^ for the 
short description of such a rule. Of course, we also 
take into account the near correlation in transverse 
momentum by limiting the maximum transverse mo- 
mentum difference for quarks and antiquarks as they 
combine into hadrons. The flavor SU(3) symmetry 
with strangeness suppression in the yields of initially 
produced hadrons is fulfilled in the model 



3 Pt spectra of the constituent quarks 
at hadronization 

It is known that the measured hadron p T spec- 
tra in relativistic heavy ion collisions exhibit a two- 
component behavior. The spectra take an exponen- 
tial form at low p T and a power-law form at high 
p T . Based on parton- hadron duality, the transverse 
momentum spectra of constituent quarks just before 
hadronization should also have the same property. In 
principle, the quarks just before hadronization come 
from two parts, i.e. the thermal quarks from the 
hot medium produced in collisions and the minijet 
quarks from initial hard collisions. The final hadrons 
are the total contribution of the two parts. But, just 
as shown in the second figure in Ref. IS, the mmi- 
jet quarks dominate the large transverse momenta 
where thermal quarks take a very small proportion, 
and thermal quarks dominate the low transverse mo- 
mentum region where the minijet quarks have a small 
contribution. Therefore, we can neglect the two small 
contributions, similar to the treatment in Ref.^J, and 
adopt a piecewise function to describe approximately 
the transverse momentum distribution of constituent 
quarks: 



— =9{p -p T )N th f th (p T ) 

2n Pt^Pt 

+0(pT-Po)N jet f ]et (p T ), (1) 

where 9(x) is the step function, N th is the number 
of thermal quarks and iV jet is the number of mini- 
jet quarks. p is the transition point from thermal 
distribution to power-law distribution, which is de- 
termined by the spectra continuity. In fact, the inter- 
action between the thermal quarks and the minijet 
quarks leads to a smooth spectrum around p , and 
we neglect this effect in this paper. 

The hot and dense quark matter produced in 
nucleus-nucleus collisions at RHIC energies shows a 
significant collective character^— '^J. The p T spec- 
tra of thermal quarks at hadronization can be de- 
scribed by a thermal phenomenological model incor- 
porating the transverse flow of thermal medium^. 
The quarks and antiquarks transversely boost with a 
flow velocity profile /3 r (r) as a function of transverse 
radial position r. (3 r (r) is parameterized by the sur- 
face velocity (3 S : (3 r (r) = (3 s £ n , where £ = r/R max , 
and R max is the thermal source maximum radius 
(0 < £ < 1). The transverse flow of thermal medium 
can be equivalently described by a superposition of a 
set of thermal sources, each boosted with transverse 
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rapidity p = tanh 1 [3 r ^ : 
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where A is the normalization constant in the region 
Pt G [0,Po]- Iq and K 1 are the modified Bessel func- 
tions. m T = \/pt 2 + m 2 is the transverse mass of the 
constituent quark. T is the hadronization tempera- 
ture. The average transverse velocity can be written 
as 



(fir) 



(3) 



With fixed hadronization temperature T = 170 MeV 
and parameter n = 0.5, the average transverse veloc- 
ity (P r ) is able to characterize the transverse collective 
flow of the hot and dense quark matter. 

The quarks and antiquarks with high transverse 
momenta are mainly from the minijets created in ini- 
tial hard collisions among nucleons. Here the so- 
called minijet quarks are those just before hadroniza- 
tion. They are the parton remnants after the revo- 
lution of the initial hard partons by gluon radiation 
and split, and are different from those in the frag- 
mentation model. The p T spectra of minijet quarks 
at hadronization can be parameterized as follows: 



/W(Pt) : 



dn 



jet 



2ir p T e?p T 



= B 1 



pr 
Po 



(4) 



where B is the normalization constant in the region 
p T € (po,oo). 

There are four independent parameters in Eq. ([1]) 
to determine the transverse momentum distributions 
of quarks: N th , N jet , {(3 T ) and a. Here, we extract the 
values of these parameters for the light and strange 
quark p T spectra at midrapidity from the data of it 
and i^°lii>i3 ; respectively. In our quark combination 
model, removing the resonance decay contributions 
from the measured ir° and K ° transverse momentum 
distributions, we get the initially produced ir° and K° s 
transverse momentum spectra. The values of param- 
eters are inversely extracted from these initial spec- 
tra. We obtain four groups of results corresponding 
to the centrality bins 0-10%, 20-40%, 40-60%, and 
60 — 80%. They are shown in Fig. 1. The lines in the 
figure are the parameterized results, from which we 
can get the quark p T spectra in any collision central- 
ity. In Fig. 2, we also show the quark p T spectra in 
four collision centralities mentioned- above. 
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Fig. 1. The values of the parameters Nth (a), 
Njet (b), {p r ) (c) and a (d) for the pt spec- 
tra of light quarks (filled circles) and strange 
quarks(open circles) at midrapidity in four dif- 
ferent centrality bins. The corresponding cen- 
trality bins are - 10%, 20 - 40%, 40 - 60%, 
and 60 — 80%. The solid and dashed lines are 
the parameterized results for light and strange 
quarks respectively. 
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Fig. 2. The pt spectra of the light quarks 
(a) and strange quarks (b) at midrapidity in 
Au + Au collisions at ^/snn — 200 GeV. 

One can see from Fig. 1 that the centrality de- 
pendence of the parameter N th is different from that 
of Nj et . As we know, the thermal quarks which dom- 
inate the low p T region carry most of the collision en- 
ergy. The available energy used to produce hadrons in 
different centrality is proportional to the number of 
nucleon participants N part at the fixed collision en- 
ergy. Therefore, the number of thermal quarks is 
approximately proportional to N part . The minijet 
quarks with high p T are mainly from the hard-jet 
created in initial hard collisions among nucleons, and 
their quantity is mainly determined by the number 
of binary collisions iV coH , which is obviously different 
from N pa rf Generally, the higher the collision cen- 
trality is, the bigger the bulk volume for hot medium, 
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the stronger the transverse collective flow (3. The 
bigger the bulk volume for hot medium is, the more 
energy loss for minijet quarks, thus the steeper the 
minijet-quark spectra and the bigger the parameter 
a. It is further observed that the strange quarks are 
different from light quarks not only in the quantity 
(due to strangeness suppression) but also in the mo- 
mentum distribution. It is easy to understand that 
the spectrum of strange quarks at high p T , due to its 
heavier effective mass, is steeper than that of light 
quarks, i.e. a (s) > a (u ' d) . In the low p T range, how- 
ever, the spectrum of strange quarks is flatter than 
the light quarks because > /3^"' d '. By analyzing 
the data of multi-strange hadrons <fr, 5 and f2, Ref. 
^ also draws the same conclusion. Furthermore, the 
similar property is also obtained in longitudinal ori- 
entation ^3. As we know, the expansion evolution 
of the plasma in the partonic phase is also a process 
of obtaining the effective mass for partons. Due to 
the heavier effective mass, the strange quarks may 
undergo a stronger hydrodynamic expansion in the 
partonic phase than the light quarks. 

As is known to all, though the nucleus-nucleus col- 
lisions at top RHIC energy exhibit a high degree of 
transparency, there are still a few net-quarks which 
stopped in the midrapidity region. The p T spec- 
tra of both 7T° and K° can not reflect the informa- 
tion of net quarks^ . We obtain the number of net 
quarks at midrapidity in different centrality bins by 
fitting the rapidity densities of net-proton^. Note 
that the ratios of tt~ /tt + and p/p measured by STAR 
and PHENIX Collaborations reveal weak dependence 
of centrality and transverse momentum^' ^ . The 
transverse momentum distribution of net quarks is 
taken to be the same as that of the newborn light 
quarks in the model. 

With the input, we can give the transverse mo- 
mentum distributions of various hadrons in different 
collision centralities and the nuclear modification fac- 
tors Rcp for these hadrons. Just as mentioned in the 
above section, we consider the decay contributions 
from all available decay channels of all resonances by 
using the decay program of PYTHIA 6.1 ^. There- 
fore, we can directly compare our calculated results 
with the experimental data. 

4 Results and discussions 

4.1 The p T spectra of hadrons in Au+Au col- 
lisions 

We firstly calculate the transverse momentum 
spectra of 7r and p (p) in different centrality bins in 
Au+Au collisions at y/s NN = 200 GeV. The results 
are shown in Fig. 2. Here, the pion spectra are cor- 



rected to remove the feed-down contributions from 
K° and A(A). The p T spectra of hadrons in low p T 
region are specially shown in the inserted plots. One 
can see that the calculated results agree well with 
the data from the STAR Collaboration. The strange 
hadron production can better reflect the property of 
the hot and dense quark matter produced in colli- 
sions. We also compute the transverse momentum 
distributions of strange hadrons K ± , A(A) in differ- 
ent centrality bins. The results are shown in Fig. 3 
and compared with the data. 
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Fig. 3. The transverse momentum spectra of 
n ± and p(p) at midrapidity in different cen- 



trality bins in Au+Au collisions at ^snn ~ 
200GeV. The data are taken from STAR Col- 
laboration [iil ]. 
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Fig. 4. The transverse momentum spectra of 
strange hadrons in different centrality bins in 
Au+Au collisions at ^/s^T = 200 GeV. The 
data are taken from PHENIX and STAR Col- 
laborations [HEl]. 
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As shown above, the good agreement between the 
model predictions and the experimental data confirms 
the validity of our model from the central to the pe- 
ripheral collisions. 

4.2 The nuclear modification factors R C p for 
hadrons in Au+Au collisions 

The nuclear modification factors Rcp can reflect 
more precisely the variation of hadron p T spectra in 
different centrality bins, ft is quantified as^: 

, _ [(d 2 N/(27T PT d PT dy))/N bm ] Ce " tral 
CP[PT) l(d2N/(2irp T dp T dy))/N Hn } p e" Phe r al ■ 

(5) 

In Fig. 4, we give the computed results of Rcp for 
7T++7T- andp+p (0-I0%/60-80%), K± and A + A 
(0-5%/60-80%), and compare them with the experi- 
mental data. 
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Fig. 5. The nuclear modification factors Rcp 
for identified hadrons in Au+Au collisions at 
y/SNN = 200 GeV. The data are taken from 
STAR Collaboration HEUEIHi- 

The hadrons in low p T region are mainly from the 
thermal quark combination. As shown in the first 
figure, the transverse collective flow of the hot and 
dense quark matter, denoted by (/3 r ), decreases with 
the falling centrality. This leads to the p T spectra 
of thermal quarks becoming softer in peripheral colli- 
sions, and then results in an increasing trend of Rcp 
for hadrons in low p T region. The hadrons with in- 
termediate p T are mainly produced by the combina- 
tion of thermal quarks with minijet quarks. As we 
all know, with the decreasing collision centrality, the 
volume size of the hot and dense quark matter be- 
comes small and the energy loss of minijet quarks 
correspondingly becomes small as they traverse the 
hot and dense medium. It leads to the p T spectra 
of minijet quarks becoming harder in peripheral col- 
lisions, which is embodied in the decreasing value of 
parameter a shown in the first figure with the falling 
centrality. With the increasing of p T , the hadron R C p 
goes over from the rise caused by the combination of 
thermal quarks to fall caused by the combination of 
minijet quarks. One can see that the quark combi- 
nation model well describes the behavior of Rcp for 
final hadrons in the whole p T region considering the 



variation of the transverse collective flow of the hot 
and dense quark matter. 

In addition, the data show that the Rcp for 
baryons clearly exhibits less suppression compared 
with that of mesons in intermediate transverse mo- 
mentum region. This type dependence of Rcp, which 
is dependent upon the number of constituent quarks 
rather than hadronic mass, has been qualitatively dis- 
cussed earlier in Refs. ^'^J as an experimental sup- 
port of the quark combination picture. Our results 
further manifest this baryon-meson difference of R C p 
in intermediate p T region. 

The experimental data show that the baryon- 
meson difference of Rcp disappears at higher p T . Us- 
ing the fragmentation mechanism, the Duke group in 
Ref. S has explained this common degree of suppres- 
sion for both baryons and mesons at high p T . Our 
results from the quark combination can also describe 
this behavior of Rcp at high p T . These two different 
hadronization mechanisms produce similar results. It 
suggests that the disappearance of baryon-meson dif- 
ference of Rcp at high p T is not caused by hadroniza- 
tion mechanism. The hard scatterings which take 
place near the surface of the collisions produce the 
back-to-back dijets. One-side jets escape almost with- 
out energy loss, while the away-side jets lose signifi- 
cant energy as they traverse the hot and dense matter. 
The hadrons with high p T in all collision centralities 
are mostly from these jets without energy loss. The 
transverse momentum distributions of these minijet 
quarks with high p T in central collisions, except the 
quantity, are almost the same with those in peripheral 
collisions. Therefore, no matter what the hadroniza- 
tion mechanism is, there exists a similar suppression 
of the Rcp for baryons and mesons at high p T . 

It is also observed that the calculated Rcp de- 
viates from the data to a certain degree. The rea- 
son may be that some effects, such as the produc- 
tion of excited-state hadrons and final-state rescat- 
tering, are not considered currently in the model. 
As we know, a small quantity of the excited-state 
hadrons are produced in relativistic heavy ion colli- 
sions, and the yields and momentum distributions of 
final hadrons are influenced by the decay contribu- 
tion of these excited-state hadrons to a certain ex- 
tent. As the decay branch ratios of many excited- 
state hadrons are incompletely measured, the contri- 
butions of excited-state hadrons are neglected in the 
current model. On the other hand, the perfect quark- 
number scaling of hadron elliptic flow v 2 suggests that 
the influence of final-state rescattering on hadron dis- 
tribution is finite and small^, so we also neglect it 
in the work. These two effects would affect the fine 
observable Rcp- 
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4.3 The p T spectra and nuclear modification 
factors R C p for hadrons in Cu+Cu colli- 
sions 

Now, we apply the centrality(participants) depen- 
dence of parameters in Au+Au collisions which de- 
termine the p T spectra of quarks in different centrali- 
zes, to the relative smaller Cu+Cu collision system at 
y/s NN — 200 GeV. In Fig. 5, we show the calculation 
results for the transverse momentum spectra of ir° , K° s 
and hyperon A in different centrality bins and Rcp 
for these hadrons in Cu+Cu collisions at y/s NN = 200 
GeV. The p T spectra of hadrons in low p T region are 
specially shown in the inserted plots. We find that the 
results agree well with the data in the low transverse 
momentum region in all centrality bins. This implies 
that the transverse collective flow of the hot and dense 
quark matter, i.e. (j3 r ), is only a function olN part and 



irrelevant to the collision system at the same collision 
energy. The good agreement of our results with the ir° 
data at high transverse momenta also suggests that 
the energy loss of minijet quarks in Cu+Cu collisions 
is almost the same with that in Au+Au collisions with 
the same participants N part at ^/s NN — 200 GeV. It 
is consistent with the recent measurement of STAR 
Collaboration^^. The above results suggest that 
the hot and dense quark matter produced in Au+Au 
and Cu+Cu collisions at the same N part and collision 
energy has similar strong-interacting character. Of 
course, even at the same N part and collision energy, 
the initial spatial eccentricity of the overlap collision 
geometry in Au+Au collisions is obviously different 
from that in Cu+Cu collision systems. The differ- 
ence is clearly reflected by some important observa- 
tions, e.g. the elliptic flow of final hadrons and the 
global polarization of hyperon^'— 
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Fig. 6. The transverse momentum spectra of identified hadrons in different centrality bins (a) and their nuclear 
modification factors Rcp (b) in Cu+Cu collisions at a/sjvjv = 200 GeV. The data are taken from PHENIX 
and STAR Collaborations [IlED]. 



5 Summary 

Using the quark combination model, we study sys- 
tematically the centrality dependence of the trans- 
verse momentum distributions for the identified 
hadrons at midrapidity in Au+Au and Cu+Cu col- 
lisions at y/s NN — 200 GeV. The centrality depen- 
dence of the parameters for quark p T spectra is ex- 
tracted from the data of ir° and K° in Au+Au col- 
lisions. We calculate the transverse momentum dis- 
tributions of 71-*, p(p), and A(A) in five central- 
ity bins. The good agreement between our results 
and the data indicates that the quark combination 
hadronization mechanism is applicable to all collision 



centralities. The nuclear modification factors Rcp for 
TT++TT- and p+p (0-10%/60-80%), K± and A+A(0- 
5%/60-80%) are calculated and compared with the 
data. The quark combination model well describes 
the behavior of Rcp for final hadrons in the whole 
p T region considering the decrease of the transverse 
collective flow of the hot and dense quark matter from 
the central collisions to the peripheral collisions. The 
disappearance of the baryon-meson difference of Rcp 
at higher p T is derived from the same transverse mo- 
mentum distribution of minijet quarks between the 
central and the peripheral collisions rather than the 
hadronization mechanism. Furthermore, we apply 
the N part dependence of parameters to the relative 
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smaller Cu+Cu collision system at the same colli- 
sion energy. We calculate the transverse momentum 
spectra of ir°, K°, and A at midrapidity in different 
centrality bins and Rcp of these hadrons in Cu+Cu 
collisions at y/s NN = 200 GeV. The results agree well 
with the data in low transverse momentum region in 
all centrality bins. It suggests that the transverse col- 
lective flow of the hot and dense quark matter is only 
the function of N part and independent of collision sys- 



tem. The calculated ir° spectrum at high p T is also 
in good agreement with the data. These results sug- 
gest that the hot and dense quark matter produced in 
Au+Au and Cu+Cu collisions at the same N part and 
collision energy has similar strong-interacting charac- 
ter. 

We are grateful to Wang Q., Liang Z. T., Yao T. 
and Han W. for helpful discussions. 
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